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Abstract 

The ground state energy and the singlet-triplet energy gap of the antifer- 
romagnetic Heisenberg model on a ladder is investigated using a mean field 
theory and the density matrix renormalization group. Spin wave theory shows 
that the corrections to the local magnetization are infinite. This indicates that 
no long range order occurs in this system. A flux-phase state is used to cal- 
culate the energy gap as a function of the transverse coupling, J±, in the 
ladder. It is found that the gap is linear in J± for J± S> 1 and goes to zero 

for J_|_ — ► 0. The mean field theory agrees well with the numerical results. 
PACS number(s): 75.30.Ee, 75.30.Kz, 75.10.Jm 
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I. INTRODUCTION 



The antiferromagnetic (AF) Heisenberg ladders (two coupled spin chains) are relevant to 
the understanding of how the physics evolves from the purely one-dimensional (ID) systems 
to two-dimensional (2D). Also, the 5=1/2 Heisenberg model on the ladder can model 
the magnetic properties of systems such as the vanadyl pyrophosphate!, (VO)2P2C>7. The 
calculation of the thermodynamic properties of this material will be addressed in a forth 
coming work!. The relevance of the transverse coupling J± for ground state properties is 
examined in the present work. 

From the theoretical point of view, two main reasons are behind the increasing interest 
in the Heisenberg ladder. The Haldane conjecture! that the energy gap of the elementary 
excitations of chains depends on whether the spin is integer or half-integer. It is well knownB 
that the spin one-half ID Heisenberg model is gapless. However it is not clear how the gap 
behaves when a transverse coupling, J± is turned on between the two chains of the Heisenberg 
ladder. In Ref. 5 the authors reported that a finite transverse coupling J± c 0.4 is required 
to get a finite gap. Latter on, Barnes et aM concluded that the critical value vanishes 
J±c = 0- The second reason is the discovery of high critical temperature superconductors 
(HT C S). The interest for these systems is due to the belief that the 2D Heisenberg model 
describes the AF interactions in the undoped copper-oxygen planes of HT C S. The ladder 
problem can help understand the crossover from ID to 2D systems (this question is related 
to the stability of the Luttinger liquid). 

Strong and have also recently used this type of model to investigate the com- 

petition between magnetic ordering and the Kondo effect in heavy fermion systems. Since 
the work of Doniachl, it is widely believed that spin excitations of heavy fermions can be 
modeled by such a model. The major weakness is that the Kondo effect is lost; only spin 
competition effect is shown. 

The problem of a plane of coupled spin 1/2 chains was studied by Azzouz! who found 
within a mean field approximation that once J± is non-zero long range AF order appears 
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in the system. The gap remains zero because of the broken symmetry due to long range 
order. For the ladder, no broken symmetry is expected to occur because of the ID nature 
of the system. Indeed, the finite size in the transverse direction will bring different physics 
than the two dimensional Heisenberg model. We believe that the gap starts to be non-zero 
for any finite J± due to the finite size in the transverse direction^. In the limit of large J± 
the system is equivalent to a weakly coupled singlets and the gap is given in leading order 
by the singlet-triplet energy separation. The first term of the gap is linear in J±. When J± 
becomes of the same order as the parallel coupling the situation becomes more difficult to 
analyze. 

In this paper, the ladder problem is investigated using a mean-field approach and exact 
diagonalization. The energy gap is found to be non-zero for any finite J±. Comparison 
with the exact diagonalization which is based on the density matrix renormalization group 
(DMRG)§ is reported. The agreement between these approaches is very good. 

II. MEAN FIELD TREATMENT 

The Heisenberg model on the ladder is denoted as follows 

h = j ± ]T Si ■ s,- + j x; St ■ s j (i) 

(*iJ>JL 

where the sums run over first nearest neighbors along the chains and (i,j)± perpen- 

dicular to the chains. J and J± are AF coupling constants. In the following J is set equal 
to unity and periodic boundary conditions are imposed. Simple limits of this model can be 
analyzed. The first obvious one is realized for Jj_ ^> 1 as mentioned in the Introduction. 
In this case, one gets weakly coupled (by J = 1) singlets and the first excited state has 
an energy gap behaving as J± in leading order. This excitation is obtained when the state 
of a single pair of spins changes from singlet to triplet. The second limit, less obvious, is 
J_i_ = 0. The two chains can be treated separately. It is known from the exact results of des 
Cloiseauxi that the ID Heisenberg model has an energy spectrum of the form 
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e(Jfe) = -| sin A; | (2) 

which shows a zero gap. The intermediate regime (J± ~ 1) is quite interesting and is the 
most complicated one. The dependence of the energy gap, hereafter denoted E g (J±), on the 
transverse coupling is investigated here using the same mean-field theory as in Ref. 7. 



A. Review of the flux-phase and Neel-flux-phase states 



The 2D generalization of Wigner- Jordan transformation of Ref. 7 can be easily imple- 
mented in the case of the ladder. One gets, following the notation of Fig. 1 

i-l 



S { x = Ci t i exp 



in ( n>t,i + ne,2 
e=o 



(3) 



for the chain 1 and 



& 2 = d , 2 exp 



i i— 1 

/7r i E n ^,i +J2 n t,2 



(4) 



for the chain 2. The indices i run along the chains. The Hamiltonian is now written in this 
fermion representation. One finds the following spinless interacting fermion Hamiltonian : 



H 



-J 



E 



+ c i)2 e 



i+5,2 



J 



^ 7T~ E c i,l°i,2 



+J E (n i)J --l/2)(n i+ i li -l/2) + J ± X;Ki-l/2)K2-l/2) (5) 

i,j=l;2,<5 i 

The phases $ are as follows 

= -7rrii_i j2 

(6) 

$ i)i+ i(2) = vrni+1,1 

$i,i-l(2) = -7T7li,i 

and 5 refers to the first nearest neighbors of a given site. The mean field solutions stud- 
ied here are the fiux-phase0, with zero magnetization and the Neel-flux-phase with finite 
magnetization m = 2(n»j) — 1. The flux due to the XY term (the first and second terms 



in Eq. ([5])) of the Hamiltonian is taken to be 7r per plaquette on average. For the Ising 
term (the last two terms of Eq. (|5|)), one chooses (c it jc\ +S j) = \(cijc\ +S j)\e~' l6i ' i + s ^ and 
( c ij4,j+s) = \( c i,j c lj+s)\ e ~ ie ' J ' J+s(l) (J = M)- In the following, we set \{c i;j 4 +SJ )\ = Q and 
I ( c hj c l,j+s) I = P- The sum over #'s around one plaquette is also taken to be 7r on average. 
The bipartite character and the different phases on each link of the system are summarized 
in Fig. 2. Despite the fact that a finite magnetization in a system like the ladder is not pos- 
sible because of its smallness, we will discuss the Neel-flux-phase and compare its physical 
meaning with the more physical flux phase state. 

B. Results and discussion 

The mean field Hamiltonian is written as follows 

H — — 5Z( C !,l e "^+1,1 + C !,1 C «-1,1 + C !,2 C «+1,2 + c l,2 et7Tc i-l,2 
J 



+ Y^( C ll C ^2 + cl,2Ci,l) 



+^ ( mn M ~ mn i+s,j + m 2 /2 \ + ^ ^(mn it i - mn it2 + m 2 /2 

i,j=l,2,S i 
i,i=l,2 V 

+J± E (Hi c ^ + Pcl^i + P 2 ) (7) 

i 

where a bipartite lattice due to AF correlations is used, Fig. 2 (the local magnetization is 
staggered: m = uia = —nis where A and B are two adjacent sites). We get 

H = YE±(k)al ± a K± (8) 

k 

in k-space where the fermionic operator a-^ is obtained from Ck through the diagonalization 
of the Hamiltonian (0). The dispersion relation is given by 

E±(k) = ± ( m 2 (l + J ± /2) 2 + (1 + 2Q) 2 sin 2 k x + (1 + 2P) 2 ( J ± /2) 2 cos 2 k y J (9) 

where J± is divided by a factor 2 since the periodic boundary conditions used in the trans- 
verse direction count J± twice. The minimization of the free energy with respect to m, Q 
and P gives a set of three self-consistent equations which become 



m = J (d 2 k/2ir)m(l + J ± /2)/E + (k), 

Q = J (d 2 k/2ir)(l + 2Q) sin 2 k x /E+(k), (10) 

P = J (d 2 k/2n)(l + 2P){JJ2) cos 2 k y /E+(k) 
at zero temperature. The integration is over the first Brillouin zone. By definition, we write 
/ d 2 k/2ir = J (dk x /2ir)(l/2) J2k y where k y can take two values: or n. One easily notes that 
m = is a solution. An interesting feature shown by such a solution is that when J± = 
the k-dependence of the dispersion relation yields 

£7±(Jfe) = ±(l + 2Q)|sinJfe|. (11) 

The ground state corresponds to the situation where the lower band is fully occupied and 
the upper band is empty. A fermion a created in the upper band produces the elemen- 
tary excitation in the system and the corresponding energy excitation is given only by the 
dispersion relation of the upper band, namely: 

e(Jfe) = (l + 2Q)|sinfc|. (12) 

One can calculate Q and finds 1 + 2Q w 1.63. This result compares well with ir/2 ~ 1.57 in 
the exact solution of Eq. @. The interesting feature is that one recovers smoothly the ID 
limit of the dispersion relation by taking m = 0. The energy gap E g (J± = 0) is then equal 
to zero. When J± is nonzero, the gap has the form 

E g (J ± )=E + (k=(0,7r)) 

= (m 2 (l + J ± /2) 2 + (1 + 2P) 2 (J ± /2) 2 ) 1/2 (13) 

which reduces to 

E,(J±) = { -^ 1 Jt (14) 

for m — 0. Eq. ([13D is obtained by calculating the difference between the ground state 
energy 

Eos = JQ 2 + J ± P 2 + (J + J ± )^ - i / |V(k) 

4 2 J 2tx 
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and the first excited state energy 



E E x = JQ 2 + J±p 2 + (J+ J±) 



2 

m 



1 r d 2 k 

2 ik^(0,7r) 27T 



E + (k) + -E + (k=(0,n)). 
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For m = the results of the numerical calculation for the set of Eqs. (|T0| ) are displayed 



in Figs. 3 and 4. The parameters Q and P show no simple dependence on J±. The 
energy gap, which is displayed in Fig. 4, has a linear behavior in J± ^> 1. This is in good 
qualitative agreement with the simple limit J± = oo. It has a more complicated dependence 
for intermediate transverse coupling because of the J^-dependence of P (Fig. 3). For small 
Jj_, E g (J±) has a simple power law form 



where the constant c = 0.76 and the exponent g = 1.15. This result compares qualitatively 
well with that of Strong and Millisi who studied this problem in the case of z-anisotropy in 
the parallel Heisenberg coupling. 

The Neel-flux-phase state has a nonzero m for J± < 1.76. The different parameters of 
this state are displayed in Fig. 5. The gap is found to go to a finite limit when J± — > 0. The 
finite magnetization in this state implies broken rotational symmetry. Gapless collective 
modes related to spin wave excitation would then exist in the gap. The spin wave theory 
goes beyond the mean field approximation. For nonzero m, the quantum fluctuations due 
to spin wave excitations would have drastic repercussions on the value of m. As in ID, 
these fluctuations destroy long range order. Indeed, for the ladder, the corrections to the 
local magnetization in the standard spin wave theory can be calculated and are found to be 
logarithmically singular 



This implies that spin wave theory is not self-consistent and no long range order can occur 
at zero temperature. So the flux-phase solution, m = 0, is more adequate to describe the AF 
correlations even if its ground state energy is slightly higher than that of the Neel-flux-phase 
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(15) 




(16) 



state as shown in Fig. 6. All the physical information the nonzero solution contains is that 
the AF correlations are more important for < J±_ < 1.76 as we see on Fig. 5 because the 
magnetization is zero for J± > 1.76, but, to our opinion, these correlations are not strong 
enough to induce long range order, (imagine that we can solve exactly this problem, then one 
would find a zero magnetization for any value of J±). The Neel-flux-phase and flux-phase 
states give the same dispersion relation for J± > 1.76. 



III. NUMERICAL INVESTIGATION WITH THE DMRG 

We have used the recently introduced density matrix renormalization group method0 
to find the ground state energy of the Heisenberg ladder. Our Fortran codes are written 
with the version of the infinite lattice method with open boundary conditions. We first 
find the ground state wavefunction of the finite 2x7 system and start the renormalization 
process by keeping 80 states in each block. Unlike the one-dimensional version suggested 
by White^ we insert only one pair of new sites in the middle of the two blocks to form 
the new superblock in each renormalization procedure so that the size of the Hilbert space 
is kept within our computer's capacity. We find this works reasonably well for the ground 
state energy shown as crosses in Fig. 6. For the singlet-triplet energy gap however, it works 
less well. Nevertheless, the gap we obtain numerically shows convincingly that it vanishes 
only as the J± goes to zero as should be apparent from the crosses of Fig. 4. The numerical 
results will be given in detail in a forth coming workll 



IV. CONCLUSION 

In conclusion, our mean-field theory (flux-phase) describes accurately the low lying ex- 
citations of the AF Heisenberg model on the ladder. The gap is found to increase smoothly 
with J±. Its behavior as a function of J± is shown in Fig. 4. The quantitative agreement 
between the analytical approach and the DMRG numerical solution is very good. 
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The accuracy of the flux-phase in the case of the ladder is a precursor of high dimension- 
ality physics since a finite flux can exit only in dimensions higher than 1. The spin energy 
excitations have a gap. Now, when charge degrees of freedom are introduced, the situation 
becomes more complicated. However, if we assume that a finite doping, 5 C , is required to 
bring the gap to zero, then one can conclude that the ID Luttinger liquid state is unstable 
for 5 < 5 C in the sense that the spin correlations decrease algebraically in ID rather than 
exponentially for finite J±. The system does not belong to the same universality class when 
J± > as that of the ID system. 
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FIGURES 
FIG. 1. The ladder system. 

FIG. 2. The flux per plaquette is equal to it on average 

FIG. 3. The parameters Q and P as a function of J± for m = 

FIG. 4. The energy gap as a function of J±. The full and dashed lines are respectively from 
our mean field treatment (m = 0) and the result of perturbation theory of Ref. 6. The +'s are the 
result of the D.M.R.G calculation. 

FIG. 5. The parameters of the nonzero magnetization solution of the mean field approximation 
(plotted as a function of Jj_) 

FIG. 6. The ground state energy as a function of J±. The full and dashed lines correspond 
respectively to zero and nonzero magnetization. The +'s are the result of the D.M.R.G calculation. 
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